INFRARED ECHOES REVEAL THE SHOCK BREAKOUT 
OF THE CAS A SUPERNOVA 



Eli Dwek^ and Richard G. Arendt ^ 
ABSTRACT 

The serendipitous discovery of infrared echoes around the Gas A supernova 
remnant the Spitzer sateUite has provided astronomers with a unique opportu- 
nity to study the properties of the echoing material and the history and nature of 
the outburst that generated these echoes. In retrospect, we find that the echoes 
are also clearly visible as infrared "hot spots" in IRAS images of the region. The 
spectra of the echoes are distinct from that of the dust in the general diffuse 
interstellar medium (ISM) revealing hot silicate grains that are either stochasti- 
cally heated to temperatures in excess of ~ 150 K, or radiating at an equilibrium 
temperature of this value. We show that the maximum luminosity that can be 
generated by the radioactive decay of ^^Ni is not capable of producing such spec- 
tra, and could therefore not have given rise to the echoes. Instead, we find that 
the echoes must have been generated by an intense and short burst of EUV-UV 
radiation associated with the breakout of the shock through the surface of the 
exploding star. The inferred luminosity of the burst depends on the amount of 
attenuation in the intervening medium to the clouds, and we derive a burst lu- 
minosity of ~ 1.5 X 10^^ Lq for an assumed H-column density of 1.5 x 10^^ cm^^. 
The average H-column density of the IR emitting region in the echoing clouds is 
about 5 X 10^'' cm~^. Derivation of their density requires knowledge of the width 
of the echo that is sweeping though the ISM, which in turn is determined by the 
duration of the burst. A burst time of ~ 1 d gives a cloud density of ~ 200 cm~'^, 
typical of dense IR cirrus. 
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ISM: interstellar dust - Infrared: general 
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INTRODUCTION 



Short-lived luminous sources can produce echoes of their outburst in the ISM. These 
echoes can be manifested as line emission from the gas, or reflected and thermally-reradiated 
light from the dust in the ISM. They can be used to probe the morphology of the interstellar 
medium (ISM) through which they are expanding and to reconstruct the historical record of 
the temporal behavior of the light source. 



^Visible light e choes have been d etected around Nova Persei 19 01 (lKapteynlll902l : lBode et al, 

2004h . V838 Mon JSond et al.lboosh . and SN 1987A JCrottJl988l l Fluorescent echoes arising 
from the equatorial ring around SN1987A have been used to probe the early hours of the burst 



from the SN, the radius of the ring, and the d istance to the LMC flFransson fc Lundqvist 



19891 : IPanagia et al.lll99ll : iDwek fc Feltenlll992f ). The SuperMACHO project has discovered 
three additional sets of light echoes in the LMC which can be traced back to known super- 
nova remnants (S NRs), and were used to derive ages of 610 and 410 yr for two of the SNRs 
fiRest et al.l[2005[ ). 



One of the most unexpected discoveries made by Spitzer satellite has been the serendip- 
itous detection of infrared (IR) light echoes in the vicinity of the Cas A SNR. The possibility 
of detecting I R light echoes from an cient SNe in general, and Cas A in particular, has been 
considered by Ivan den Berghl (119651 ) more than 4 decades ago. Cas A is the remnant of a 
supernova that exploded ~ 320 yr ago in the Perseus Arm of the Galaxy, 3.5 kpc away. The 
echoes were first (and best) detected in 24 /xm MIPS scans of Cas A (Fig. [1]). Subsequent 
observations revealed that the echoes are also visible in all IRAC channels and ground-based 
K band observation s. The echoes vary on time scales of < 6 months, and extend at least 
1.5° from the SN R jHines et al.l 120041 : iKrause et al.l[2005l : iKim et al.ll2008[ ). More recently, 
Kim et al.l (120081 ) have used the 24 fim maps of the echoes over a three year interval to probe 
the structure of the ISM around Cas A, revealing its complex filamentary structure. 

In this paper we combine early spectral observations with more recent data, that have 
revealed a wider extent of the echoes, to provide a detailed analysis of the nature of the 
radiative burst that gave rise to the echoes, and the density of the interstellar clouds giving 
rise to the echoes. 

We first present a brief review of the discovery of the echoes around Cas A, and derive 
the IR spectrum of select echoing clouds spanning a wide range of angular distances from 
the SNR (§2). We also provide a brief mathematical review of the geometry of an echo, 
and derive the distances of the echoing clouds from the center of the explosion. In §3 we 
present the model input parameters: the burst delay time and its spectral characteristics, 
and the interstellar dust model used in the calculations. A stochastic heating model is used 
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to calculate the thermal emission from the dust exposed to the burst of radiation, and in 
§4 we describe the dependence of the calculated IR fluxes on the burst characteristics. The 
model fits are then used to derive the properties of the burst and the echoing clouds (§5). 
The results of our paper are briefly summarized in §6. 



2. THE INFRARED ECHOES AROUND CAS A 

2.1. The Discovery of the Echoes 

The initial MIPS 24 /xm observations (a 0?2 x 0?3 scan map) of Cas A serendipitously 
revealed "chains of IR knots" extending up to ~ 12' from the center of the SNR (Hines et 
al 2004). The team followed up with ground-based K band (2.2 /xm) observations which 
revealed extremely high proper motions of the features. This motivated a second epoch of 
MIPS 24 /im observation that confirmed that much of the small scale structure outside of the 
Cas A SNR was changing on time scales of several months. A subset of the echoes seemed 
to represent a bipolar jet, nearly oriented in the plane of the sky, leading to the suggestion 
that at least these echoes m ay have been produ ced by a recent flaring of the Cas A neutron 



star around 1952±2.5 A.D. (iKrause et al.l 120051 ). i.e. about 50 yr ago. 



Subsequent MIPS 24 /im observations have been repeated at roughly 6 month intervals, 
and have been expanded to include at first a ~ 1 square degree field around the SNR, and 
later a ~ 9 square degree field. Examination of these data from the Spitzer archive reveal 
that echoes are not restricted to the bipolar structure originally detected. They are found at 
all position angles with respect to the SNR and extend at least 1?8 from the SNR (limited by 



the siz e of the map). The neutron star burst scenario was dismissed in passing by lKim et al. 



( I2OO8I ). Here we will provide additional physical evidence against such scenario. 



Follow-up observations of the echoes also included low and high resolution IRS obser- 
vations of a bright echo in a filament located 0?8 from the center of the SNR. Additional 
low resolution spectra were collected for outlying regions in the vicinity of the SNR. Two of 
the closest regions to the SNR, located within the NE jet of ejecta and the SW counter jet, 
have emission line s and the 22 nm silicate feature that are characteristic of ejecta material 



( lEnnis et al.ll2006l : iRho et al.ll2008l ). Spectra of other regions exhibited a fiat continuum and 
coincide with highly variable features, and are therefore assumed to be representative echo 
spectra. 



^Archival MIPS images were obtained from Spitzer programs PID = 718, 30571 ~ G. Rieke, and PID = 
231, 233, 20381 - O. Krause. 
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The discovery of the echoes with the Spitzer led us to examine if the echoes could also 
be identified (in hindsight) in the map of the Cas A region that was obtained in 1983 with 
the Infrared Astronomical Satellite (IRAS). Figure [2] (left panel) shows a false color image at 
12, 25, and 60 /im of the reprocesses IRAS map of Cas A and its environment. The emission 
is dominated by cool interstellar dust grains, peaking at wavelengths of ~ 100 — 140 /xm 
(red regions in the image). However, there are several regions with excess 25 /zm emission 
(including Cas A itself) from hot dust, which thus appear green. Comparison to the 24 /xm 
MIPS data (right panel), smoothed to /iJ^S" resolution, clearly shows that the 25 /im excess 
features are echoes which have moved between 1983 and 2003. The fact that many echoes 
are in similar locations indicates that many of the clouds producing the echoes are structures 
that extend for at least 20 lyr along the line of sight (cf. Kim et al. 2008). IRAS and Spitzer 
have observed different portions of these large-scale structures. 



2.2. The IR Spectrum of Select Echoes 

For the purpose of this work we have concentrated on the analysis of six IR knots 
that were identified as echoes and have been the target of IRS observations. The selected 
echoes, listed in Table 1 and shown in Figure [T|, sample a range of angular distances from 
the nominal center of the explosion. Also included in the table is a seventh 24 /xm echo, 
without a spectrum, as a representative of the most distant echoing cloud. The SL+LL IRS 
spectra of the echoes were collected from the Spitzer archive^] Because there are often several 
sources or extended emission within the slits, we have reprocessed the basic calibrated data 
(BCD) using the SSC tools IRS_CLEANMASK and SPICE. The off-target subslit positions 
were used to subtract background emission for the corresponding on-target position. In a 
few cases where the off-target subslit position crossed the SNR, a clean subslit position from 
an alternate target was used. Both SL and LL data were processed, although in most cases 
the SL integration times were too short to detect the echoes, which are visible in all IRAC 
bands. Therefore the analysis presented here relies only on the LL spectra. Figure [3] presents 
the ~ 14 — 40 /im spectra of the echoes listed in Table 1. 

All echoes within a 1000" of Cas A exhibit a steep rise in their spectrum between 
~ 14 — 20 /xm, not seen in the IR spectra of dust iri the diffuse ISM that is heated by the 



general interstellar radiation field (jPraine fc Lil 120071 ). This spectral characteristic provides 



important clues on the nature of the emitting dust, strongly constraining the viability of 



^Archival IRS spectra were obtained from Spitzer programs PID = 3310 - L. Rudnick, and PID = 20381 
- O. Krause. 
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different echo models. Echo 6 is brighter than the echoes closer to the SNR, and is part of 
a large extended structure that can be seen in the older IRAS data. This region appears 
to contain a non-variable emission component which would simply be the same cloud(s) 
illuminated by the diffuse interstellar radiation field. We have therefore subtracted a scaled 
ISM spectrum from the Echo 6 spectrum such that its 20 - 35 fim slope is within the range 
of the slopes of the other echo spectra. This non- variable component accounted for 20 - 50% 
of the 20 - 35 /xm emission. 



2.3. The Location of the Echoes 

We assume that the echoes are the reradiated thermal emission from dust heated by the 
UV- visual output from the SN, located at point S (see Figure Hj). The emission from dust 
located at point C will be detected by an observer at point O with a delay time t given by 
ct = r + X — d, where r and x are, respectively, the distances of C from the source and the 
observer, and d = OS is the distance of the observer to the source. The locus of all points 
with equal delay time is an ellipsoid with the source and the observer at the focal points. 

Using the law of cosines: = + rf^ — 2xdcos{a), we can express the distance r in 
terms of the delay time t and source distance d as: 

^_d{d + ct){l-cos{a)) + {cty/2 

^ ~ d{l- cos(a)) + ct ^ ' 

where a is the angular distance of the echoing dust from the source. 

Figures [S^-c are three different depictions of the interrelations between the delay time, 
and the angular and physical separations of the echoing cloud s from the source assuming a 
source at the distance of Cas A, corresponding to d = 3.4 kpc ( Reed et"al1ll995l ). Figure [5^ 



shows the echo distance from the source, r, as a function of the delay time t for the separation 
angles a corresponding to the echoes listed in Table 1. All distances are bound by the 
lower limit ct/2 which is the minimum distance corresponding to dust located at point A 
immediately behind the source (see Figure H]). The vertical lines reprsents the maximum 
delay time, corresponding to t = 320 yr, the age of Cas A, and a delay time of 50 yr. 
Figure Eb depicts the echo distance from the source as a function of angular separation. The 
vertical dashed lines correspond to the angular separations of the echoes listed in Table 1. 
The filled squares and open diamonds correspond to distances calculated for delay times of 
ct = 320 and 50 yr, respectively. Figure (St depicts the physical location of the echoes with 
respect to the source, as projected on a plane that includes the source, located at the origin 
{x,y) = (0,0), and the observer (off scale), located at {x,y) = {d,0). The vertical line going 
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through the origin represents the plane of the sky going through the source. The two eUipses 
shown in the figure are the loci of all points with equal delay times. The filled squares and 
open diamonds represent the echoes listed in Table 1 assuming delay times of 320 and 50 yrs, 
respectively. Echo 7 with the largest separation angle from the source is almost on the plane 
of the sky when t = 320 yr, but at a large distance in front of the source when t = 50 yr. 

The figures illustrate that when delay times are short, a given range of angular distances 
will correspond to a wide range of echo-source distances, whereas long delay times will 
translate ino a significantly narrower range of echo-source distances for the same range of 
angular distances. When t = 50 yr, the range of distances, r, spanned by the echoes in 
Table 1 is between 27 and 1172 lyr. This range is narrowed to only 160 and 359 lyr when 
t = 320 yr (see Table 2). 



3. MODEL INPUT PARAMETERS 

3.1. The Echo Delay Time 

The most obvious explanation for the echoes is that they were generated by the radiative 
energy released following the cor e collapse of t he pr ogenitor star of the Cas A SN, which 



occurred about 320 yr ago. Since iKrause et al.l (120051 ) suggested that a subset of the echoes 



were produced by a recent fiaring of the neutron star that happened around 1950, we will 
also examine if such outburst could have produced the other echoes as well. 

For each of the two scenarios {t = 50 or 320 yr) we determine the incident fiux required 
to produce each echo spectrum. Since the two scenarios place the echoing dust clouds at 
different distances, these fluxes will translate into different burst luminosities. To check the 
viability of any physical scenario we require that all echoes are generated by a single burst 
with a well-deflned intrinsic luminosity. 

It is possible to construct scenarios involving multiple bursts from a neutron star over 
the past 320 years. However, such scenarios have too many free parameters to be usefully 
constrained by the observations. Furthermore, the construction of such scenarios seems 
unnecessary given the consistent results that are ultimately obtained with the 320 yr old SN 
echo scenario. 
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3.2. The Burst Spectrum and Duration 



Because neither the Cas A SN nor any subsequent outburst of the neutron star has been 
observed, we do not know the spectrum of the radiation that gave rise to the echoes. We 
will therefore explore three possible "burst" spectra: (1) an EUV spectrum characterized 
by a temperature of T;, = 5 x 10^ K; (2) a UV spectrum characterized by a temperature of 
Tb = 5 X 10^ K; and (3) an optical spectrum characterized by a temperature of = 6000 K. 

The radiative burst giving rise to the IR echo could be generated by the breakout of 
the SN shock wave through the stellar surface, in which case the burst is powered by the 
kinetic energy of the explosion. A flash of EUV-UV shock breakout radiation has been 
predi cted by hydrodynamical simulations of SN explosions (IKlein &: Chevalierl 1 19781 : iFalk 
19781 ). however, because it lasts for only a few days, it has only recently been directly 
observed in a SN. Its existence was first indirectly in ferred in SN1987A from the detection 



of na rrow UV and optical emission lines from the SN ( IFransson et al.lll989l : ISonneborn et al. 
19971 ). The lines arise from an equatorial ring, located about 0.7 lyr from the SN and 



produced by r nass loss from th e prog enitor star that was flash-ionized by the EUV burst. 
More recently, ISoderberg et al.l (120081 ) . using the Swift X-ray Telescope, detected an X-ray 
burst from SN 2008D which they attributed to the shock breakout. 

Alternatively, the IR echo could be powered by the energy released in the decay of ra- 
dioactive elements that have formed in the SN explosion. In addition to the optical emission, 
a signiflcant UV component has bee n observed in the light curves of all major types SNe 
with the Swift UV-optical Telescope ( iBrown et al.ll2008l ). 



The IR spectrum of the echo depends on the spectral shape and intensity of the flux 
incident on the dust. The luminosities of the burst are free parameters that will be derived 
from model flts to the IR echoes. The duration of the burst is constrained by energy con- 
siderations and only affects the thickness and overall brightness of the expanding light echo. 
The burst duration is a necessary parameter in converting the derived column density into 
a volume density of the radiating dust. 



3.3. Interstellar Dust Model 



We used the BARE-GR-S model of IZubko et al.l (120041 ) to characterize the echoing 
interstellar dust. The model consists of a population of bare silicate and graphite grains 
and PAH molecules (see Table 5 in their paper). In this model, the dust-to-gas mass ratio is 
0.0062, and the abundance of its various dust constituents is consistent with solar abundances 
constraints. The size distribution of the different dust species are depicted in Figure 19 of 
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their paper, and the parameters of their functional forms are hsted in their Table 7. 



4. MODELING THE ECHO SPECTRA 
4.1. The IR Emission from Stochastically Heated Dust 

The observed IR flux, F^{X), from an extended source of radiating dust is given by: 



0.2 



f{a) k{v, a) ^y(A, a) da 



0.1 



(2) 



where Vt is the angular size of the dust cloud, is the mass column density of the radiating 
dust, /(a) is the grain size distribution normalized to unity over the {01,02} size interval, 
and k(z/, a) is the mass absorption coefficient of the dust. Implicit in eq. ([2]) is a sum over 
dust composition. The term fi;(z/, a) ^^(A, a) in the equation is the spectral response of the 
dust to an incident radiation field, with Qy given by an integral over dust temperatures T^: 



g,{X,a)= P{T,)nB,{Ta,X)dT, (3) 

where P{Td) (K~^) is the temperature probability distribution of the stochastically heated 
dust, and B^{Td, A) is the Planck function. For dust radiating at the equilibrium temprature 
rpeg^ P(Trf) collapses to the delta-function S(T^'^), and Qi^{\,a) = 'iTB^(T^'^,a) 



4.2. Echo Spectra For Different Burst Models 

To model the IR emission from the echoing clouds we will expose a population of in- 
terstellar dust to a burst of radiation emanating from the SN. The burst spectrum will be 
characterized by a blackbody at a fixed temperature Th. A significant fraction of its soft 
X-rays and EUV photons will be absorbed in the intervening medium to the clouds. The 
optical depth, t(z/), will be parametrized by the intervening ISM column density, N^- The 
flux incident on the dust grains is then given by: 

where Fq is a fixed fiux, calculated for a nominal luminosity Lq = 1 x 10^^ Lq and a nominal 
distance of tq = 160 lyr between the radiating dust and the source, and ^ is a fiux scaling 



7iB,{Tb, A) 



-r(A) 



(4) 
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parameter determined by the fit to the observed spectrum of the IR echo. Numerically, Fq 
is given by: 

Fo = = 1.33 X 10^ erg s"^ cm'^ . (5) 
47rro 

Since the H-column density to the echoing clouds is unknown, we will calculate the 
burst's attenuation for three different H-column densities given by Nh = 1.5 x 10^^, 1.5 x 10^^, 
and 1.5 X 10^° cm~^. For a typical echo distance of 160 lyr these column densities correspond 
to average ISM densities of uh = 0.01, 0.1 and 1.0 cm~^, respectively. 

Our burst model calculations were performed for seven different burst spectra character- 
ized by different burst temperatures and intervening H-column densities. Burst models are 
designated by their spectra: EUV for = 5 x 10^ K; UV for T;, = 5 x 10^ K; and OPT for 
Tj, = 6000 K, and by a suffix of 18, 19, or 20, corresponding to the assumed SN-echo column 
density. For example, model EUV19 corresponds to a burst with a temperature of 5 x 10^ K 
and a column density of 1.5 x 10^^ cm^^ to the echoing clouds. Model UV20 corresponds 
to a burst with a temperature of 5 x 10^ K and an ISM column density of 1.5 x lO^'' cm~^. 
We neglected any attenuation for the OPT burst. This model lacks any significant ionizing 
radiation, and therefore the extinction was negligible even for the largest column density. 

We calculated echo spectra for each burst spectrum for a grid of fluxes, characterized by 
different values of ^ ranging from 10~^ to 10. For example, we calculated 16 IR spectra for 
burst model EUV18 by varying ^ from 0.001 to 1.0 with 14 intermediate values. The only 
free model parameter is the normalization factor, which is equal to the dust mass column 
density, J^d, of the model spectrum that produced the lowest value of x^- 

Figure E] compares the fluxes from the different burst scenarios outlined above as they are 
incident on an echoing cloud at a distance of 160 lyr from the SN. For illustrative purposes, 
the fluxes from the EUV and UV bursts were calculated for values of = 1.0, and the flux 
from the optical burst was calculated for ^ = 0.01, reflecting its expected lower luminosity. 
For sake of comparison, the figu re also includes the flux of the general interstellar radiation 



field [ISRF; iMathis et al.l fll983f )] 



Figures [71 [8] and [9] depict the results for the OPT, the EUV and the UV bursts, respec- 
tively. The top panels in the figures show the calculated IR spectra of dust that is exposed 
to a range of illuminating fluxes, characterized by different values of ^. To illustrate the de- 
pendence of the spectral signature of the dust on the flux and spectral shape of the incident 
radiation we normalized all fluxes to unity at 20 /im. The middle panels depict the best fits 
to the echo spectra by the incident fluxes. Different echoes required exposure to different 
fluxes to fit their IR spectra. The bottom panels of the figures show the contributions of the 
different dust components to the total flux from Echo 3. 
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The figures show that the radiation from the optical burst cannot reproduce the observed 
characteristic rise in the echo spectra from 14 to 20 fim. In contrast, both the EUV and UV 
models reproduce the echo spectra fairly well, independent of the density of the intervening 
medium. 

The decomposition of the fits to Echo 3 into the various dust components illustrates 
that the rapid rise in the echo spectra between ~ 14 — 20 /im, is caused by the 18 fim silicate 
feature. Two conditions have to be met to produce the observed rise in the echoes' spectra: 
(1) the sihcate grains have to be heated to temperatures in excess of ~ 150 K; and (2) the 
silicate emission has to dominate the emission from graphite grains which have a smooth 
spectrum over the ~ 10 — 40 fim wavelength interval. The echoes generated by the optical 
burst fails to meet both conditions. Silicate grains have a low opacity to optical photons and 
therefore remain at relatively low temperatures. The top panel of Figure [TD] shows that the 
equilibrium silicate temperature is only about 83 K for the OPT burst. Stochastic heating 
is limited to the smallest grains with sizes < 70 A. In contrast, graphite grains have a much 
larger opacity at visible wavelengths and attain much higher equilibrium temperatures of 
~ 175 K. Consequently, the IR echo spectrum produced by the optical burst is dominated 
by graphite emission. 

The EUV and UV bursts meet both requirements for producing the rise in the echo 
spectra. As shown in Figure [TOl they both possess sufficiently hard photons to raise the 
equilibrium silicate temperature to ~ 150 K. Compared to the optical burst, a larger fraction 
of the small grains are fluctuating to temperatures beyond ~ 100 K. Furthermore, the opacity 
of silicate grains rises dramatically in the EUV, and becomes comparable to that of the 
graphite grains. Consequently, as illustrated in the bottom row of Figs. [8] and [H the total 
IR spectrum is dominated by sihcate emission. 



4.3. A Fit to the Average Echo Spectrum 

It is of interest to examine to what extent our model fits to the ^14-40 nm echo spec tra 



also fits the average spectral energy distribution (SED) obtained by lKrause et al.l (120051 ) in 
the 2.2, 3.6, 4.5, 5.8, and 24 /im bands. This SED, shown in Figure [HI was obtained by 
averaging over 59 position in the northern lobe region of the preliminary MIPS image around 
the remnant. The bold line in the figure shows the model fit to the IRS spectrum of echo 3 
(orange line). The X symbols represent this model's fiux densities after integration over the 
IRAC and MIPS 24 fim bandpasses. These fiux densities were normalized to fit the average 
echo fiux density at 24 /im. The fit to the average echo spectrum required the abundance 
of PAHs to be increased by a factor of 2 over their general ISM abundance. The increased 
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PAH abundance did not affect the fit to the spectrum of echo 3 since PAHs contribute 
relatively little e mission at wavelen gths > 14 fim. Also shown in the figure is the average 



cirrus spectrum (IZubko et al.ll2004l ). 



The figure shows that if the echoes are only observed in the broad IRAC and MIPS 24 fim 
bands, their colors are very similar to those interstellar cirrus dust. Consequently, the true 
temperature of the echoes is totally lost with this limited photometric data. In contrast, the 
IRAS 25-to-60 fj,m color ratio was clearly able to distinguish between the echoing clouds and 
the interstellar cirrus (see Fig. [2]), confirming the unique spectral signature of the echoes. 
The temperature of the echoing dust is much higher than that of the dust in cirrus clouds, 
requiring an incident flux that is significantly harder than the general ISRF. 

The excess 2.2 fim emission is caused by reflection, and may therefore sample a totally 
different portion of the SN light curve, powered by the radioactivitis in the ejecta. The 
intensity of the reflected echo could therefore provide valuable information of the Cas A 
progenitor. The column density of the scattering dust could therefore be significantly higher 
than that of the thermally emitting dust, providing a geometry-dependent observational test 
of this hypothesis. 



5. MODEL RESULTS 

5.1. Burst Scenario and Derived Burst Properties 

The incident flux that provides the best fit to a given echo spectra can be converted to 
burst luminosity by: 

^6=(^) ^fLo (6) 

where r is the actual distance of the particular echo from the SN, Lq = 1 x 10^^ Lq, and 
is the value of ^ that produced the best fitting incident flux. 

Each echo requires a particular incident flux to produce its IR spectrum. Given this 
flux, the value of depends on the assumed column density of the intervening ISM [see eq. 
(jlj)]. A higher column density will require larger values of ^p, corresponding to larger burst 
luminosities. 

Figures [T2l and [T3] depict the burst luminosities required to produce the echo spectra for 
the different burst models and delay time scenarios. We assume that all echoes were produced 
by a single burst. A viable echo scenario should therefore give the same burst luminosity 
for all echoes. The results show that a scenario with the 50 yr delay time produces a large 
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spread in burst luminosities, compared to the 320 yr delay scenario. The figures shows that 
the dispersion in the derived burst luminosities for the latter scenario is only about 50% 
of the mean value (see also Table 3), compared to a factor of two in the former. This is a 
direct consequence of the fact that the echo spectra are generally very similar, thus requiring 
exposure to the same incident fiux of radiation. Narrowing the range of burst luminosities for 
the 50 yr delay model will require more distant echoes to have systematically lower column 
density, a no n-physical a nd co ntrived solution to this problem. We therefore support the 



conclusion of iKim et al.l (120081 ) that the echoes were not generated by a hypothetical burst 
from the neutron star occurring about 50 yr ago. In principle, all echoes could have been 
produced by repeated bursts occurring over a period of 320 yrs. However, a more likely 
scenario is that the all echoes were produced by a single burst associated with the Cas A SN 
event. 

Table 2 lists the model fits for each of the echoes assuming a 320 yr delay time. For each 
model we present the incident fiux required to produce the best fit to the echo spectrum as 
characterized by the value of C,f, the burst luminosity. Lb, the column density of the echoing 
cloud, Nh, and the of the fit. For each echo, the value of was normalized to unity 
for the best fitting model. From the values of the of the fit, we see that the OPT burst 
provides a significantly worse fit to the echo spectra than the EUV and UV bursts. The 
EUV and UV bursts are generally of similar quality, although for echoes 1, 2, and 3, UV 
models are distinctly better. The models with a lower H-column density usually provided 
slightly better fits than the models with Nh = 1.5 x 10^° cm~^. 



The lack of any evidence for H-recombination lines from the gas around Cas A (IKrause et al 



20051 ) suggests that the burst ran out of ionizing photons before it reached the echoing clouds. 
Inspection of Figure [6] shows that these conditions are easily met when the intervening H- 
column density is sufficiently high with A^^^^ > 10^^ cm^^. 

Table 3 lists the burst luminosity averaged over all echoes for each burst model. It 
clearly shows the correlation between the derived burst luminosity and the column density 
of the intervening ISM. Larger column densities require higher burst luminosities to produce 
the same fiux needed to produce the observed echo spectrum. Excluding model results for 
the highest column density, we get that the average EUV-UV burst luminosity is = 
1.5 X IO^Lq. 
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5.2. The EUV-UV Burst Must be Powered by Shock Energy 



In principle, the EUV-UV luminosity could be either powered by the explosion energy 
or by the radioactive decay energy in the ejecta. A simple analysi s, based on the best 
estimate of the properties of the Cas A progenitor (lYoung et al.ll2006l ). firmly rules out the 
latter possi bility Based o n exp losion calculations, ejecta velocity and mass constraints, and 
^^Ti yields, lYoung et al.l (120061 ) estimate the mass of ^^Ni to be between 0.06 and 0.2 M©. 
Taking the average decay energy to be 1.73 MeV and a decay lifetime of 8.8 d, the maximum 
luminosity that can be powered by the decay of ^"^Ti is between 1.2 and 4.1 x 10^ L©. This 
luminosity is about two orders of magnitude smaller than that required to produce the IR 
echoes. 

The luminosities required to produce the echoes are easily obtained by tapping into 
the kinetic energy of the explosion. Taking an average explosion energy of 2 x 10^^ erg 



(lYoung et al.l 120061 ) suggests that this energy must be released during a period of 40 d 
to produce an average luminosity of 1.5 x 10^^ L©. Several hydrodynamic model s for the 
explosion of SN1987 A were constructed to fit the earli est UV and opti cal light curve (IWoosley 
19881 : lEnsman fc Burrow]ll992l : [Blinnikov et allboooh . The models of iBhnnikov et all kood l 
show that the burst had a peak blackbody temperature of ~ 10^ K, cooling on a timescale 
of days to a blackbody temperature of ~ 6000 K, the H-recombination temperature in the 
expanding photosphere. Burst durations are therefore typically < days. Taking a burst 
decay time of one day, suggests that only 2.5% of the kinetic energy of the explosion need 
to be thermalized to produce the required EUV-UV luminosity. These simple arguments 
conclusively show that the echoes have revealed the breakout of the SN shock through the 
stellar surface. The combined constraints on the energy of the explosion, the required burst 
luminosity, and hardness of its spec trum can provide useful constraints for modeling the 
structure of the progenitor of Cas A (IMatzner fc McKed 119991 ). 



5.3. Derived Cloud Properties 

From the fits of the calculated IR spectra to the observations we can derive the column 
density, of the emitting dust. We take fl to be the width of the IRS LL slit squared, 
fl = (10.6")^. The corresponding H-column density is given by: N h = -Md/im H Zdu) , 



where Zan = 0.0087 is the dust-to-H mass ratio of the BARE-SG-S IZubko et all (120041 ) 
interstellar dust model, and tuh is the mass of the H-atom. The column densities derived for 
each echo and burst model are listed in Table 2. Echo 6 is the brightest, and consequently 
has the highest H-column density. Table 4 lists the column density of each echo, averaged 
over all EUV and UV burst models. The average column density for Echoes 1-5 is about 
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4.5 X 10^^ cm ^, and higher by a factor of ~ 5 for Echo 6. 



All echoes within 50' of Cas A are located behind the SN. The IR emission arises from 
a narrow region of thickness cAth, where Atf, is the burst duration, along the observer's 
line of sight that is illuminated by the burst of radiation. The width of the radiation front 
is determined by the burst duration At;,. The hydrogen number density, nn, of the echoing 
cloud can be derived from the dust mass column density, A^^, by: 

Md \ ^ ^_l 385 



cm 



-3 



(7) 



where the numerical value was derived for a nominal echo column density of A^^ = 1 x 
10^8 cm-2. 

The burst duration should be taken as the length of time d uring which the burst s pec- 
trum is dominated by EUV-UV emission. From the models of iBlinnikov et al.l (120001 ). we 
take this time to be about 1 d, giving an average cloud density of uh ~ 385 cm~^. This 
i s a plausible densit y for the echoing clouds, corresponding to those of typica l cirrus clouds 



( Wolfire et al.l 120031 ). Cas A lies on the far side of the Perseus arm [e.g. iBieging et al. 



19911)]. so there are re latively few (if any) denser molecular clouds behind the remnant 



( lUngerechts et al.ll2000l ) to reflect the radiative burst from the SN. 



6. SUMMARY 

The Spitzer satellite discovered a series of IR echoes around Cas A that are caused by 
the delayed arrival of thermal emission from dust heated by the radiative output from the 
SN. Echo spectra are distinct from the IR emission from the diffuse ISM, exhibiting a sharp 
rise at wavelengths from ~ 14 to 20 /xm. In hindsight, the echoes can also be found in the 
IRAS images from 1983. 

We assume that the echoes are generated by thermal emission from interstellar dust 
heated by a single burst of radiation released after the SN event. We calculated the fluxes 
incident on the dust required to generate the observed echo spectra for different bursts 
characterized by blackbody spectra with temperatures of 5 x 10^ (EUV burst), 5 x 10^ (UV 
burst), and 6000 K (optical burst). 

Our models show that the rise in the echo spectra is due to the 18 /im silicate feature. 
The silicate grains radiate mostly at the equilibrium temperature of ~ 120 K and dominate 
the echo spectra. The optical burst is incapable of heating the silicates to sufficiently high 
temperature required to reproduce the spectral characteristics of the echoes, and is therefore 
ruled out as the source of the echoes. 
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The echoes are instead generated by the intense EUV-UV burst of radiation generated 
by the shock breakout. We show that the echo spectra can be reproduced by bursts with a 
wide range of temperatures from ~ 5 x 10^ to ~ 5 x 10^ K. The burst luminosity required 
to generate the echoes depends on the assumed H-column density between the SN and the 
echoing cloud. We find the average burst luminosity to be ~ 1.5 x 10^^ Lq for an average 
intervening H-column density of ~ 1.5 x 10^^ cm~^. 

The average H-column density of the echoing clouds is ~ 5 x 10^'' cm~^, which for a 
burst duration time of ~ 1 d gives an average cloud density of ~ 200 cm^"^. This density is 
typical of dense cirrus clouds, and consistent with the location of Cas A behind the Perseus 
arm and the paucity of dense molecular material behind the remnant 

The luminosity required to produce the echoes cannot be generated by the radioactive 
decay of the 0.06-0.2 Mq of ^^Ni that has formed in the explosion, and must therefore reflect 
the fraction of the kinetic energy of the explosion that was converted to EUV-UV radiation 
as the SN shock broke out through the stellar surface. The Cas A echoes represent the first 
indirect "view" of a shock breakout via the thermal dust emission from echoing clouds. 

ED acknowledges the support of NASA's LTSA03-0000-065. We thank the referee for 
his/her critical comments. This work is based on observations made with the Spitzer Space 
Telescope, which is operated by the Jet Propulsion Laboratory, California Institute of Tech- 
nology, under contract with NASA. 
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Table 1. Observed Properties of Echoing Clouds 



Target 


Location (RA: DEC) 


a cr 


'"so (lyr)^ 


^"320 (lyr)^ 


Echo 1 


23h23'38.57": 


+58d52'54.4" 


259" 


27.0 


160.3 


Echo 2 


23h22'50.58": 


+58d47'30.6" 


300" 


27.6 


160.4 


Echo 3 


23h22'41.55": 


+58d48'37.0" 


359" 


28.7 


160.6 


Echo 4 


23h22'50.98": 


+58d41'49.4" 


508" 


32.5 


161.2 


Echo 5 


23h24'12.37": 


+59d01'12.8" 


820" 


44.5 


163.1 


Echo 6 


23h21'39.96": 


+59d34'26.9" 


2860" 


258 


198 


Echo 7 


23h26'38.62": 


+57d01'19.4" 


6630" 


1172 


359 



^Angular distances were calculated from the optical expansion 
center located at (RA : DEC) = (23h23'27.77": +58d48m49.4s) 
( jThorstensen et al.ll200ll ) 

^Echo distances from the SN for a delay time of t = 50 yr. 

^Echo distances from the SN for a delay time of t = 320 yr. 
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Table 2. Derived Model Properties^ 



Echo # 




1 


2 


3 


4 


5 


6 


OPT: 




0.6 


0.6 


0.8 


0.2 


0.3 


0.6 




Lb 


0.60 


0.60 


0.81 


0.20 


0.31 


0.92 




Nh 


2.7 


1.6 


1.9 


1.3 


3.0 


9.7 






2.20 


2.46 


4.62 


1.81 


1.72 


1.31 


EUV18: 


^F 


0.2 


0.2 


0.2 


0.04 


0.08 


0.1 




Lb 


0.20 


0.20 


0.20 


0.041 


0.083 


0.15 




Nh 


4.7 


2.8 


4.2 


4.4 


7.1 


32.0 






1.67 


1.66 


3.98 


1.00 


1.03 


1.13 


EUV19: 




0.2 


0.2 


0.4 


0.06 


0.1 


0.2 




Lb 


0.20 


0.20 


0.40 


0.06 


0.10 


0.31 




Nh 


7.7 


4.6 


3.5 


4.9 


9.5 


26.9 




x' 


2.40 


2.94 


3.21 


1.05 


1.40 


1.00 


EUV20: 


Cf 


2.0 


2.0 


2.0 


0.4 


0.8 


1.0 




Lb 


2.0 


2.0 


2.0 


0.41 


0.83 


1.53 




Nh 


3.5 


2.1 


3.1 


3.2 


5.2 


23.8 




x' 


1.93 


1.74 


3.16 


1.03 


1.00 


1.10 


UV18: 


Cf 


0.1 


0.1 


0.2 


0.02 


0.06 


0.1 




Lb 


0.10 


0.10 


0.20 


0.020 


0.062 


0.15 




Nh 


5.8 


3.5 


2.8 


5.4 


5.9 


20.6 




x' 


1.00 


1.52 


1.00 


1.64 


1.10 


1.09 


UV19: 


Cf 


0.2 


0.2 


0.2 


0.04 


0.08 


0.1 




Lb 


0.20 


0.20 


0.20 


0.041 


0.083 


0.15 




Nh 


3.6 


2.2 


3.2 


3.2 


5.3 


24.0 




x' 


1.14 


1.00 


1.45 


1.28 


1.07 


1.29 


UV20: 


Cf 


0.4 


0.2 


0.4 


0.06 


0.1 


0.2 




Lb 


0.40 


0.20 


0.40 


0.061 


0.10 


0.31 




Nh 


3.0 


3.4 


2.7 


3.5 


6.9 


20.0 




x' 


1.50 


2.23 


1.17 


1.41 


1.61 


1.22 



is the value of ^ [see eq. ^] that gives the best fitting 
incident flux [eq. Burst luminosities are in units of 10^^ Lq, 
and echo column densities are in units of 10" cm . For each 
echo, is renormalized to 1.00 for the best fitting model. 



-20- 



Table 3. Average Burst Luminosities 



Burst Model 


<Lb> 


a 




(1012 Lq) 


(1012 L^) 


OPT 


0.57 


0.28 


EUV18 


0.15 


0.070 


EUV19 


0.21 


0.13 


EUV20 


1.47 


0.70 


UV18 


0.11 


0.064 


UV19 


0.15 


0.07 


UV20 


0.25 


0.15 



< Lb > is the value of Lb averaged over 
all echoing clouds for each burst model, 
and a is the r.m.s. dispersion of the value. 



Table 4. Average Echo Column Densities^ 



Echo # 


<Nh> 


a 




(1017 cm-2) 


(1017 cm-2) 


1 


4.7 


1.8 


2 


3.1 


0.9 


3 


3.3 


0.6 


4 


4.1 


1.0 


5 


6.7 


1.6 


6 


24.5 


4.4 



i< Nh > is the value of Nh averaged 
over all EUV-UV burst models for each 
echo, and a is the r.m.s. dispersion of the 
value. 
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Fig. 1. — A map of the region around Cas A showing the location of all transient light signals identified 
as echoes (red squares). The echo enclosed by a white rectangle is also shown as an inset in the upper right 
corner of the figure. The inset shows the position of the IRS SL and LL slits. 
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Fig. 2. — Left panel: This false color image displays the region around Cas A as observed by IRAS in 
1983. Reprocessed IRAS data (IRIS) at 60, 25, and 12 fim are depicted as red, green, and blue respectively. 
The extended structures with excess 25 fim emission (green in color) are now know to be light echoes as they 
have significantly changed in Spitzer MIPS images in 2006. Cas A is the bright green source at the center of 
the image. Right panel: This image shows IRAS 25 fim observations from 1983 in green, and Spitzer MIPS 
24 /im observations from 2006 in magenta. The MIPS data have been smoothed to approximate the IRAS 
spatial resolution. Cas A is the bright compact source at the center of the map. The light echoes are visible 
as extended structures, with magenta or green tints depending on when they were present. Nonvariable 
structure appears a uniform gray shade. Most of the 1983 (green) echoes seen here are also evident via their 
mid-IR colors as revealed the the previous figure. 
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Fig. 4. — The geometry of an echo. The locus of all points with equal delay time is described by an 
ellipsoid with the source "S" and the observer "O" at the focal points. The lengths AS = ct/2, and BS = ct, 
where t is the difference between the observing and the emission time of a given photon. The point "C" is 
an arbitrary point on the ellipse sustending an angular distance a from the source, located at distances r 
and X from the source and observer, respectively. The length OS = d, the distance of the source from the 
observer. 
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Fig. 5. — Three different depictions of the interrelations between the delay time, and the angular and 
physical separations of the echoing clouds listed in Table 1 from the source. The relation between the 
different quantities is given by eq. ^ for d — 3.4 kpc, the distance to Cas A. Top: Echo distance to the 
source versus the delay time for different angular separations. Middle: echo distances to the source as a 
function of angular distance for delay times of 50 and 320 yr. Bottom: Physical location of the echoes with 
respect to the source, as projected on a plane that includes the source, located at (0,0), and the observer 
(off scale), located at (d, 0). All figures are described in more detail in §2.3 of the text. 
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Fig. 6. — The flux of the different bursts considered in this paper (§4.2): the EUV burst characterized by 
a blackbody temperature of Tb = 5 x 10^ K (top panel), and the UV burst with = 5 x 10'' K (bottom 
panel). Also shown in the figure is the flux from the optical burst with Th = 6000 K. For comparison, we 
also included the flux of the ISRF. The intrinsic luminosities of the EUV and UV burst were taken to be 
10^^ Lq. The luminosity of the optical burst was taken to be 10^° i©. All fluxes were calculated for a 
distance of 160 lyr The unattenuated burst is shown as a violet line, and the solid, dotted, and dashed line 
represent the flux attenuated en route to the echoing cloud by the ISM with column densities of 1.5 x 10^*, 
1.5 X 10^^, and 1.5 x 10^° cm^^, respectively. Detailed burst parameters are described in the text. 
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Fig. 7. — Model fits for the optical burst (OPT). Top row: Calculated IR spectra of an echo exposed to 
different fluxes F* characterized the product ^Lq and calculated for = 160 lyr (see eq. |4]) with ^ — 0.1, 
0.3, 0.5, 0.7, 0.9, and 1.0. All fluxes were calculated for rg = 160 lyr (see eq. S]) and are normalized to unity 
at A = 20 /im. Middle row: Best fit models to the IR echoes. Bottom row: Decomposition of the best 
fitting IR spectrum of Echo 3 into its emission components from silicate and graphite grains and PAHs. The 
figure illustrates the failure of the OPT burst to fit the echo spectra. A detailed discussion of the figure is 
in the text. 
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Fig. 8. — Model fits for the EUV burst propagating tlirough a densities of 0.01 (left column), 0.10 (middle 
column), and 1.0 cm^^ (right column). Top row: Calculated IR spectra of an echo exposed to different 
fluxes characterized by values of ^ = 0.001, 0.008, 0.06, 0.10, 0.40, and 1.0 (left panel); ^ = 0.001, 0.006, 
0.02, 0.08, 0.40, and 1.0; and £_= 0.1, 0.4, 0.8, 2.0, 6.0, and 10. All fluxes were calculated for tq = 160 lyr 
(see eq. H|) and are normalized to unity at A = 20 ^m. Middle row: Best fit models to the IR echoes. 
Bottom row: Decomposition of the best fitting IR spectrum of Echo 3 into its emission components from 
silicate and graphite grains and PAHs. A detailed discussion of the figure is in the text. 
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Fig. 9. — The same as Figure |8] for the UV burst. Calculated IR spectra of an echo exposed to different 
fluxes characterized by values of ^ = O.Of , 0.04, 0.08, 0.20, 0.60, and f .0 (left and middle panel); and ^= 
0.04, 0.08, 0.20, 0.60, f .0, and 4.0. All fluxes were calculated for rp = i60 lyr (see eq. |4]) and are normalized 
to unity at A 20 /im. Middle row: Best fit models to the IR echoes. Bottom row: Decomposition of 
the best fitting IR spectrum of Echo 3 into its emission components from silicate and graphite grains and 
PAHs. A detailed discussion of the figure is in the text. 
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Fig. 10. — The temperature probability distribution P{Td) (K^^) for 0.0010 fim (dashed hue); 0.0030 
(dashed-dotted Une), and 0.010 /im (sohd Une) sihcate grains exposed to the optical burst (top panel), the 
EUV18 burst model (middle panel), and the UV18 burst model (bottom panel) with fluxes calculated for ^ 
= 0.6, 0.2, and 0.1, respectively [see eq. [4])]. Larger grains radiate at the equilibrium dust temperature. Only 
the EUV and UV burst are capable of heating a significant fraction of the silicate grains to the temperatures 
(> 100 K) that are needed to produce the 18 fim feature. A more detailed discussion of the figure is in the 
text. 



- 31 - 



1000.00 I 



100.00 



10.00 - 



■-5 




10 100 

A(yu.m) 



1000 



Fig. 11. — The spectral energy distribution of echoes in the K-band (2.2 /xm), IRAC (3.6, 4.5, 5.8, and 
8.5 fim) and MIPS (24 /ini) bands, averaged over 59 positions in the northern lobe region of the MIPS image 
Krause et all |2005l ) is represented by filled green circles. The X symbols in the figure represent our model fit 
(bold solid line) to the spectrum of echo 3 (red line), convolved with the IRAC and MIPS 24 /im bands, and 
normalized to the 24 /im MIPS datum point. The fit to the average echo spectrum required the abundance 
of PAHs to be increased by a factor of 2 over their general ISM abundance. Also shown in the figure is 



the average cirrus spectrum (jZubko et al.ll2004l ). Observed only in the IRAC and MIPS 24 /im bands, the 
average echo spectrum is very similar to that of interstellar cirrus. The excess 2.2 /im emission is caused by 
direct reflection a different component of the SN lightcurve. 
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Fig. 12. — The EUV burst luminosities required to generate the IR echoes from Cas A [see eq. ([6])] for 
different values of the density of the ISM around Cas A: uh = 0.01 cm^'^ (top panel); nn = 0.1 cm^"^ 
(middle panel); and nn = 1.0 cm~'^ (bottom panel). Filled squares are values calculated for a delay time of 
320 yr, whereas the open diamonds are calculated for a delay time of 50 yr. A viable model should give a 
single value for Lh for all echoing clouds. The scenario in which the delay time is 360 yr is more physical, 
since it produces a significantly narrower range of burst luminosities compared to the 50 yr model. The 
spread in burst luminosities in t=320 yr scenario may reflect the variations in the burst attenuation en route 
to the echoing cloud. 
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Fig. 13. — Same as Figure [12] for the UV burst. 
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